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O . ABSTRACT 

i We calculate the abundance of dark-matter concentrations that are sufRciently over- 

' dense to produce a detectable weak-gravitational-lensing signal. Most of these overden- 

Qh, sities are virialized halos containing identifiable X-ray and/or optical clusters. How- 

ever, a significant fraction are nonvirialized, cluster-mass, overdensities still in the 
process of gravitational collapse — these should produce significantly weaker or no X- 
ray emission. Our predicted abundance of such dark clusters are consistent with the 
abundance implied by the Erben et al. (2000) detection of an apparent dark lens. Weak 
I lensing by these nonvirialized objects will need to be considered when determining cos- 

^ . mological parameters with the lens abundance in future weak-lensing surveys. Such 

weak lenses should also help shed light on the process of cluster formation. 



< 

(N 
(N 



O 

m 
o 

1 INTRODUCTION 

O ' 



Key viTords: galaxy clusters — weak gravitational lensing — cosmology 



Weak gravitational lensing due to the deep gravitational potential of a galaxy cluster gives rise to a detectable weak distortion 
of the images of background galaxies. This weak shear has now been detected around roughly 30 clusters and been used to 
map the total dark-matter mass in the clusters as well as the dark-matter distributions within the clusters (see Bartelmann 
O ' & Schneider 2001, Mellier 1999). Weak lensing also has the potential to map the mass distribution on even larger scales 
' (Miralda-Escude 1991; Blandford et al. 1991; Kaiser 1992; Bartelmann & Schneider 1992; Stebbins 1997; Kamionkowski et 
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^2 ' al. 1998). Just last year, four groups independently reported detection of cosmic shear, distortions to background galaxies 
' induced by weak gravitational lensing by mass inhomogeneities on few-Mpc scales along the line of sight (Bacon et al. 2000; 
Kaiser et al. 2000; Wittman et al. 2000; Van Waerbeke et al. 2000). It is apparent that in the future, such cosmic-shear surveys 
I J will have the sensitivity to identify galaxy clusters in the field. Since such surveys will probe the total mass directly, it could 
rS [ provide a powerful new technique for determining the cluster-halo abundance and thus the power-spectrum amplitude as and 
^ . matter density Q,m (e.g., Kruse & Schneider 1999; Reblinsky et al. 1999). 
■ ■ ■ ' In fact, one spectroscopically-confirmed cluster has already been detected via its gravitational-lensing effect on background 
galaxies (Wittman et al. 2001). More intriguing is the apparent dark lens discovered by Erben et al. (2000). This lensing signal 
corresponds to a ~ 10^*Mq mass concentration, but there is no obvious corresponding galaxy overdensity (Gray et al. 2001) 
and only faint (if any) X-ray emission. Evidence for other apparent dark lenses has been reported by Miralles et al. (2002) 
and Koopmans et al. (2000), the latter involving a detection through strong, rather than weak, lensing. 

In retrospect, the existence of such dark concentrations should not come as too much of a surprise. Galaxy clusters form 
at rare (e.g., > 3a) high-density peaks of a Gaussian primordial distribution. Thus, for every virialized cluster, there should 
be a significant number of proto-clusters (e.g., 2a — 3a peaks), mass overdensities that have not yet undergone gravitational 
collapse and virialized, but which have begun to break away from the cosmological expansion. The timescale for collapse of 
cluster-mass objects is large, and the overdensities can be very large even before they have virialized. It should thus not be 
too surprising if such objects produce a weak-lensing signal that resembles that from virialized clusters. 

These proto-clusters should contain galaxies and maybe a few groups that later merge to form the cluster. Since the X-ray 
luminosity is a very rapidly varying function of the virialized mass, the summed X-ray emission from these objects should be 
much smaller than that from a fully virialized cluster of the same mass. When we refer to these proto-clusters as "dark," we 
thus mean that they should be X-ray underluminous. Strictly speaking, the mass-to-light ratios of these clusters should be 
comparable to those for ordinary clusters. However, high-redshift clusters may be difficult to pick out in galaxy surveys, and 
these proto-clusters should have a sky density a few times smaller. Thus it would not be surprising if these dark lenses had 
no readily apparent corresponding galaxy overdensity. 

In this paper we calculate the abundance of dark and virialized lenses. To do so, we first determine the overdensity 
required to produce a detectable weak-lensing signal as a function of redshift. We consider several different density profiles 
including a homogeneous sphere, an isothermal sphere, a Navarro, Frenk, & White (Navarro et al. 1997; Navarro et al. 1996; 
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Navarro ct al. 1995) profile and a Hcrnquist (Hcrnquist 1990) profile. We then use the spherieal-top-hat-coUapse (STHC) 
model to determine the differential abundance of overdensities as a function of position along their evolutionary cycle. Using 
the aperture mass technique (Schneider 1996) we can then determine the sky density and redshift distribution of halos that 
are sufficiently overdense to produce a detectable weak-Iensing signal. 

As our results below will show, there should be rougiily one dark lens for every 5-10 virialized lenses discovered by weak 
lonsing. It is worthwhile to point out that this result is robust in that the ratio of dark to virialized lenses is not expected 
to be very sensitive to the amount of observational noise in the lensing map, i.e., observational noise will equally affect the 
detectability of both types of lenses. Therefore, although our results are obtained by assuming the only source of noise is the 
intrinsic ellipticity distribution of the source galaxies - in accordance with other such theoretical weak-lensing studies found 
in the literature - the predicted relative abundance of dark and virialized lenses will not change very much if we made a 
more exact estimate of the total noise in weak-lensing maps. It is also encouraging to note that given the sky coverage and 
average image size of weak-lensing maps to date, the number of dark lenses we would expect to have seen is of order unity 
and therefore consistent with the detection (Erben et al. 2000, Miralles et al. 2002) of one or two dark lenses. 



2 MINIMUM OVERDENSITY REQUIRED TO PRODUCE WEAK-LENSING SIGNAL 

In this Section we provide the conditions for an overdensity of mass M and radius R at redshift z to produce a detectable 
weak-lensing signal. Following the procedure of Bartelmann & Schneider (2001) (see also Schneider 1996, Seitz & Schneider 
1997, Kruse & Schneider 1999) we determine the dependence of a lensing system's signal-to-noise ratio on that system's 
overdensity and redshift. 

In a weak-lensing map, a mass overdensity causes the image of the background source galaxies to be tangentially sheared. 
Noise is introduced by both the intrinsic ellipticity of these background galaxies as well as by the presence of foreground 
galaxies in the image. To arrive at a signal-to-noise relation for a weak-lensing system, consider N galaxy images each at 
angular position 9i — {9i cos (j>i,6i sin 4>i) with tangential ellipticity et{0i) and within a lens-centered annulus that is bounded 
by angular radii ^in < 0i < Oout ■ The shear 7 is related linearly to the dimensionlcss surface mass density of the lens, which 
is the physical surface mass density S(0) divided by the critical surface mass density Ecrit- For a lens at redshift Zd and a 
source at redshift Ze, 

s„,(^.;^.) = ^^, (1) 

where Da, and Dds are the angular-diameter distances between the lens and the observer, the source galaxy and the 
observer, and the lens and the source, respectively. To account for the redshift distribution of the source galaxies, define (Seitz 
& Schneider 1997) 

, , _ lim^^^oo Scrit(Zd; «s) Scritoo (^d) 

yi(Zs;Zd) = 7- — -T = 7 r- (2) 

^crit\Zd,Zs} ZjcrityZd, Zs ) 

Then the dimensionless surface mass density is given by, 

Kie,z,) = = ^ n(e)Z{z.;zd). (3) 

Furthermore, the linear relation between the shear and surface mass density implies that they have the same dependence on 
source redshift so that 7(0, «s) = Z{zs; Zd)j{0)- For the rest of this paper any reference to k or 7 refers to k{0) and 7(6), 
respectively. Assuming the intrinsic orientation of galaxy sources is random, the expectation value of the image ellipticity is 
(Seitz & Schneider 1997; Bartelmann & Schneider 2001) 

Eie) ^ (Z)7(0), (4) 
where 

(Z) = / dzs Pz{zs)Z{zs; Zd), (5) 



and Pz{zs) is the redshift distribution of source galaxies. The function (Z) = {Z){zd) — of order unity for the redshifts 
considered — allows a source redshift distribution to be collapsed onto a single redshift Zs satisfying Z{zs) = (Z) (see Bartel- 
mann & Schneider 2001). 

Using the Map-statistics introduced by Schneider (1996), define a discretized estimator for the spatially filtered mass 
inside a circular aperture of angular radius 6, 

N 

Map = iVet(0OQ(|0i|), (6) 
n ^—^ 

where n is the number density of galaxy images and Q is a weight function that will be chosen later so as to maximize the 
signal-to-noise ratio of the estimator. Assuming the ellipticities of different images are uncorrelated the dispersion of Map can 
be obtained by squaring (6) and taking the expectation value, yielding 
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where ffe is the dispersion of the two component elhpticity. By (4) the expectation value of Map is, 

(Map) = ^f]7t(0i)Q(|ei|), (8) 

i = l 

where yt is the tangential shear. Taking the ensemble average of (8) over the probability distribution for the galaxy positions 

gives, 



(Map)c = 27r(z) / dee{-yt){e)Qie), (9) 

•'fin 

where {"ft) {9) is the moan tangential shear on a circle of angular radius 6 and the subscript 'c' stands for continuous. Similarly, 
we can take the ensemble average of the dispersion (7), to obtain 

2 P^ont 

ai = ^ j deeQ^{e). (lo) 

The ensemble-averaged signal-to-noise ratio is then, 

s (M.p)e 2{z)^reTde0{^mQ{e) ^^^^ 

By the Cauchy-Schwarz inequality the signal-to-noise ratio of the estimator is maximized if 

QWoc(7t)W. (12) 
Since 

lam = m - («>w, (13) 

(Bartelmann 1995) where (k)(6') is the dirncnsionlcss moan surface mass density on a circle of radius Q and niff) is the 
dimensionless mean surface mass density within a circle of radius d, the maximized signal-to-noise ratio becomes 



S_ 
N 




dee\ji{e)- {K){e)f. (i4) 



If instead of using a maximized weight function Q we chose one of the often used generic weight functions given in Schneider et 
al. (1998), our estimate of the signal-to- noise ratio for a given lens would be slightly smaller. In an upcoming paper (Weinberg 
& Kamionkowski 2002) we show that although using such a weight function reduces the predicted abundance of dark and 
virialized lenses somewhat, the principle result of this paper, namely that the relative abundance of dark to virialized lenses 
is 10-20%, is virtually unchanged so long as ^out ^ 3 arcmin. 

To compute the signal-to-noise ratio for a lens with a given density profile we need to determine the mean tangential shear 
of the source galaxies. Different density profiles will in general produce different shear patterns. In particular, the more cuspy 
a profile is the stronger its lensing signal. Of course this becomes more complicated when considering profiles with power-law 
breaks. For instance, although the NFW profile goes as at small radii while the isothermal sphere goes as r~^, at larger 
radii the former varies as while the latter remains at . The net effect, as wo will show, is that the NFW profile yields 
a stronger signal compared to the isothermal sphere for lenses at reasonable redshifts. That said, we consider a variety of 
profiles to account for the full range of possibilities and to study the dependence of our results on these profiles. Specifically, 
we compare the calculated abundances assuming the overdensity is a point mass, a uniform-density sphere, an isothermal 
sphere, an NFW profile and a Hernquist profile. For an object of a given mass, mean overdensity, and density profile, we can 
solve for the parameters of the given profile (e.g., the radius, the velocity dispersion, the scale radius, the scale density, etc.) 
and determine, using equation (14), whether such an object produces a sufficiently large weak-lensing signal-to-noise ratio so 
as to be detectable. Note that for an overdensity with angular radius smaller than the size of the lensing image, the shear 
pattern beyond the radius of the overdensity will be that of a point mass. Furthermore, if the angular radius is larger than the 
image size then the lensing signal is determined by just the mass Mp within the projected image radius P = OoutDd and not 
the mass outside this radius. The derivation of the signal-to-noise relation for each of these profiles is given in the Appendix. 

To produce a detectable signal, an overdensity must be large enough to yield a signal-to-noise ratio greater than some 
minimum value. For the calculations done in this paper we adopt the following fiducial values, unless stated otherwise: 
(yS/Ar)inin = 5, ^out = 5 aicmiu, the number density of galaxy images is n = 30 axcmin^^, and cTj = 0.2. The minimum 
nonlinear overdensity corresponds to a particular position along the linear-theory evolutionary cycle. In the next Section we 
discuss how we relate the minimum nonlinear overdensity to a corresponding minimum linear-thoory overdensity. This will 
enable us to apply the Press-Schechter formalism to obtain an estimate of the abundance of overdensities that produce a 
weak-lensing signal as a function of redshift. 
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3 DYNAMICS 

We use the STHC model to relate the minimum non-linear overdensity needed for a detectable weak-lensing signal at a given 
redshift to a minimum linear-theory overdensity. According to STHC the non-linear evolution of cosmic density fluctuations 
is approximated by a dynamical model in which the initial linear perturbation is an isolated, uniform sphere surrounded by 
unperturbed matter. Gravitational instability causes the initially small linear perturbation to grow and enter the nonlinear 
regime, ultimately forming a virialized object that is decoupled from the cosmological background. In order to avoid the 
collapse to infinite density predicted by the solution of STHC, we invoke a simple smoothing scheme that allows us to map a 
linear overdensity greater than the critical linear density contrast, 5c ~ 1.69, to a finite nonlinear overdensity. In what follows 
we shall consider the STHC model in a ACDM universe. Following the derivation of the relevant STHC formula, we present 
our smoothing scheme. Finally, we discuss how we distinguish "virialized" clusters from those that have not yet collapsed. 

For a flat cosmology with a cosmological constant, the change in the proper radius, r, with scale factor a for a uniform 
spherical overdensity of fixed mass M is given by (see Peebles 1984, Eke et al. 1996) 

\da) a-^+uja-^ ' ^ ' 

where a = {1 + z)^^ , /3 is a constant which is positive for overdensities and 

uj = (fio ' - 1), (16) 

where fio is the cosmological density parameter. Note that the units of r are such that {3M/4npoy^^ = 1 where po is the 
cosmological background density at ^; = 0. Separating the variables in equation (15) and integrating gives 

r.'l/2 „'l/2 

r*' = / — (17) 



(a;r'3-/3r' + l)'''' 7„ {cja'^ + l)^^"" 

Solving for the root of the numerator in equation (15) gives the turnaround radius (i.e., radius at maximum expansion), rta, 
as a function of the density parameter uj and perturbation amplitude (3. An exact solution for rta is given in Appendix A of 
Eke et al. (1996). For overdensities that axe past turnaround the left-hand side of equation (17) is integrated from zero to rta 
and added to the integral from r to rta- The evolution of the radius of an overdensity as a function of time is illustrated in 
Figure 1. Note that the cosmological constant has the effect of slowing the collapse as compared to a CDM universe. 
The non-linear overdensity is given by 

^^^NL^Ppert^ (18) 
Pb 

where Pport is the mean density of the perturbed region and pb is the background density at the given redshift. Since 
Ppert = Pa/r"^ and pb = po/o^, the non- linear overdensity becomes 

For a given non-linear overdensity of mass M at redshift z we can find the radius of the perturbation r such that (S/N) > 
{S/N)jnin- We can then solve equation (17) for /3. 

We now relate this same /? to the linear-theory perturbation amplitude. Eke et al.(1996) showed that 

^=3^(a„(2^)V3)'5» , (20) 
where ao is the scale factor today, <5q" is the linear-theory overdensity extrapolated to the present and 

I ' ^ ^ du, (21) 



(Peebles 1980). The linear-theory overdensity at redshift z is given by 

,5"°(«) = 5i'"D(a), (22) 
where D{a), the linear theory growth factor for a ACDM cosmology, is 

A{a{2c.y/^) 

^(«)=A(a„(2.)V3)- (23) 
Using equations (20), (22), and (23) we get the desired relation between the linear-theory overdensity and /3: 
un, . ^4(a(2aO^ 

^ ^"^ = ^^ ao(2a.)V3 " (2^) 

Equations (17), (19) and (24) therefore provide a map between the non-linear overdensity and the linear-theory overdensity 
at a given redshift. 

It can be shown that r — > in the limit that 5''° — > 6c, corresponding to the well known inflnite density predicted by the 
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Figure 1. The radial evolution of a density perturbation according to the STHC model. At the turnaround time, t = tta, tlie perturbation 
reaches a maximum-expansion radius and begins to collapse. As expected, in a ACDM cosmology {dashed curve) the collapse taJses 
somewhat longer than in a CDM cosmology {line-dot curve). The collapse to a singularity predicted by the solution of the STHC model 
is avoided by the smoothing scheme {solid curve) which yields a constant radius once the virialized overdensity is reached. 

solution of STHC. An actual overdensity will, of course, virialize before reaching the singular solution. To properly account 
for this we introduce the following smoothing scheme. 

Rather than assuming that an overdensity satisfies equation (19) throughout its evolution, assume it satisfies it only until 
it reaches the virialized overdensity 1 + 5^^{z). Once the perturbation reaches the virialized overdensity take its radius to 
be a constant with time so that the overdensity continues to grow only because the cosmological background density keeps 
decreasing. The non-linear overdensity is therefore given by 



where Ovir is the scale factor at virialization. Since (5''°(t2) = (5''°(ti)Z)(a2)//3(ai), the linear-theory overdensity then becomes 



Therefore, if the minimum non-linear overdensity needed to produce a detectable weak-lensing signal at redshift z is larger than 
the virialization overdensity, we evaluate Ovir using equation (25) and then compute the minimum linear-theory overdensity 
using the lower expression in equation (26). In Figure 1 we plot the radius of an overdensity as a function of time using this 
smoothing scheme. In Figure 2 we show the non-linear overdensity as a function of the linear-theory overdensity. Note that 
the value of the overdensity at virialization can be obtained by assuming r = r\ir, the virialized radius, in equation (19), 
and using the expression from Lahav et al. (1991) which gives the ratio between the turnaround radius and the virialization 
radius. For convenience we use the Kitayama & Suto (1996) approximation to 1 -I- 5^^{z), as well as their approximation to 
5c{z). We independently verified that both approximations matched the solution of the exact formalism described above. 

In summary, given the minimum non- linear overdensity needed to produce a detectable weak-lensing signal, 5™^,^, of an 
object of mass M at redshift z, we use equations (17), (25) and (26) to compute the corresponding minimum linear-theory 
overdensity, <5ji"n, needed to produce a detectable weak-lensing signal. If 5m°n < &c{z) then the object can produce a detectable 
weak gravitational lens, even though it is not yet virialized. 




(25) 




(26) 
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Figure 2. The nonlinear ovcrdensity as a function of the linear-theory ovcrdcnsity according to the STHC model. The full solution of 
the STHC model predicts collapse to an infinite overdensity as <5iin — > 1.69 {dashed curve). According to the smoothing scheme, however, 
once a mass concentration reaches the virialization overdensity l-{-S^^{z), its radius remains constant so that the overdensity increases in 
proportion to the decrease in the background density. The solid curves show the smoothing scheme solution for mass concentrations that 
reach the virialization overdensity at z = 0.1,0.2,0.5, and 1.0, from top to bottom. In an Einstein-de Sitter Universe, the virialization 
overdensity is independent of redshift and therefore all of the solid curves would be the same. 



4 ABUNDANCES 



To calculate the abundance of overdensities that produce a detectable weak-lensing signal as a function of redshift, we use 
Press- Schechter theory assuming Gaussian statistics for the initial linear-theory density field. The differential number count 
of lensing objects per steradian, per unit redshift interval is 



dV dzdQ' 



dN{3l^,^) ^ dNjC^J dV 

dzdQ, 
where 
dV c 



dzdil 



(i+^)-Av(c)^ 
Ho y/Qoil + z)^ + 1 - ^o' 



(27) 



(28) 



is the comoving-volumo clement, c is the speed of light. Ho is Hubble's constant, and Da(z) is the angular-diameter distance 
at redshift z. The total number density of weak lenses is given by 



rfAr(Cin) 
dV 



-f 

Jo 



f{M-5'X)^{M)dM, 



(29) 



where dn{M) /dM, the comoving number density of virialized objects of mass M in the interval dM, is (Press & Schechter 
1974) 



,,,, 
dM^^^ 



2 po 5c{,z) 
TV M'^ a{M,z) 



din a 



dlnM 



exp 



6,iz) 



2a2 



(30) 



In this paper we use the Viana & Liddle (1999) fits to the dispersion of the density field, a{M,z), obtained from the galgixy 
cluster X-ray temperature distribution function. The function f{M; 5m°n) is the fraction of objects, either dark or virialized. 
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Figure 3. Tlic iiicau tangential shear as a function of angular distance from the lens center for an NFW {solid curve), Hernquist {dotted 
curve) and isothermal sphere {dashed curve) density profile. The normalization is arbitrary. 



that can lens {d > (5m°n) relative to those that are virialized {5 > 5c). The probability that an object's linear overdensity is in 
the range 6i < 6 < 62 is 



P{Si <6<S2) = erf 



V2a{M, z) 



Verff 

J V V2a{M, z) 



where 'erf is the error function. Therefore, for dark lenses (i.e., those objects with 6^^^^ < & < 5c) 



/dark (M, 2) 



PjS^^i^ <5 <6c) rlin ^ c 

P{6 > Sc) ' ™ " 



(31) 



(32) 



0, otherwise, 
while for virialized lenses {6 > 5c and 6 > (5m"n)i 

'^m'in > Sc, 



/vir(M, Z) 



{ P{S > S]Sn) 

P(d>5c) ' 



(33) 



otherwise. 



For low enough masses, the minimum overdensity needed to produce a detectable weak-lcnsing signal becomes so large that 
both /dark and /vir approach zero, thereby imposing an effective weak-lensing mass threshold (see Section 5.2). Integrating 
equation (27) over redshift assuming / = /dark yields the number count of dark lenses per unit area on the sky and similarly 
for virialized lenses when / = /vir. 



5 RESULTS 

5.1 Minimum overdensity as a function of redshift 

We can now compute the sky density of weak lenses. To gain physical insight into the results as well as illustrate the 
calculational procedure discussed above we first show the redshift dependence of the minimum non-linear overdensity. As 
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Figure 4. The minimum nonlinear overdensity needed to produce a detectable weak lensing signal as a function of rcdshift for a IO^^Mq 
object with a density profile that is a uniform-density sphere {line-dot curve), a truncated isothermal sphere {dashed curve), a Hernquist 
profile {dotted curve), and an NFW profile {solid curve). An overdensity with a larger mass will displace these curves downwards. The 
thin, long-dash-dot curve is the overdensity at virialization in the STHC model. 



noted earlier, the result is sensitive to the lens density profile on account of the minimum overdensity's dependence on the 
shear. Since the shear is proportional to the surface mass density, the NFW and Hernquist profiles (for whom p oc r^^ as 
r —^ 0) have a constant shear at small radii while the isothermal sphere profile (p oc r~^) has a shear that goes as for all 
radii. This is shown in Figure 3, where we plot the radial dependence of the mean tangential shear for these different profiles. 

In Figure 4, the minimum non-linear overdensity as a function of rcdshift for a 10^'' Mq object is plotted for the various 
profiles. All the profiles show the same general trend: a minimum at a ~ 0.3 and monotonic rises at lower and higher redshifts. 
This is a consequence of the source-galaxy redshift distribution which we assume is given by a function of the form 

with /3 = 1.5 and mean redshift (zs) ~ 1.5«o = 1.2 (cf. Small et al. 1995; Brainerd et al. 1996; Cohen et al. 2000) . Since 
lenses are most effective when they lie midway between the source and the observer (i.e., the factor DdDds/Ds peaks when 
Dd — Dds), an overdensity at z ~ 0.3 is ideally positioned to lens source galaxies that are primarily located at z = {z.,} ^ 1, 
thereby accounting for the minimum in the curves. Accordingly, ovordcnsities located at lower and higher redshifts than z ^ 0.3 
are less effective at lensing so that a larger overdensity is needed to produce a detectable lens. In addition, for an overdensity 
with redshift approaching unity, there are fewer background galaxies to lens (less signal) as well as more foreground galaxies 
in the image (greater noise), further decreasing the observed lensing signal-to- noise ratio. 

Another feature to note in Figure 4 is the difference in amplitude of 1 -I- between the different profiles. Over most 
of the rcdshift range, the NFW profile requires the smallest overdensity in order to produce a detectable weak-lensing signal 
while the uniform-density sphere requires the largest. This is because the NFW profile has its mass much more centrally 
concentrated as compared to the uniform-density sphere. A source galaxy at some angular radius near the lens center will 
therefore be sheared more strongly by the former and hence produce a larger signal. A similar explanation accounts for the 
differences in amplitude of 1 -|- between the non-uniform profiles. 



5.2 The abundance of dark and virialized lenses 

In Figure 5 we show the redshift distribution (normalized to unity) of dark and virialized lenses for the NFW, Hernquist, and 
isothermal-sphere profiles. Because the minimum overdensity for the uniform-density sphere was so large, the probability of 
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Figure 5. Redshift distribution of daA lenses {left panel) and virialized lenses {right panel) for a truncated isothermal sphere {dashed 
curve), an Hernquist profile {dotted curve), and an NFW profile {solid curve). The ordinate gives the normalized probability distribution 
per unit redshift interval. 



detecting a lens with such a profile is neghgible and hence no longer considered. For all three profiles the distribution peaks 
at a « 0.5 and has a full-width at half-maximum of Az w 0.5. The distribution drops off at w 1 for two reasons: the 
minimum ovcrdcnsity is becoming increasing large since (^s) ~ 1 and the STHC dynamics predicts fewer and fewer massive, 
large overdensities at these higher redshifts. 

The sky density of dark lenses as a function of redshift for the same three profiles is shown in Figure 6. Depending on 
the density profile, we expect to find between 1 — 20 dark lenses per square degree out to 2 = 1 and virtually none at higher 
redshifts. The reason the Hernquist profile predicts a smaller dark-lens sky density compared with the isothermal sphere and 
NFW is that such a profile requires a larger overdcnsity to produce a detectable weak- lensing signal (see Figure 3). Finally, 
note that although this distribution is integrated over dark lenses of all masses, the minimum overdensity as a function of 
redshift becomes so large for M< 5 x 10^^ M© that there are virtually no dark lenses with such small masses. This point is 
illustrated in Figure 7, where we plot the weak- lensing mass distribution (i.e., the integrand of equation (29) times the mass) 
for both dark and virialized lenses at 2 = 0.5. Furthermore, since the Press-Schechter mass function falls off steeply with mass, 
there will be very few dark lenses with M> 10^^ Mq despite the lower value of the minimum overdensity at these masses. 

In Figure 5 and 6 we also show the redshift distribution and sky density of virialized lenses for the three different density 
profiles. Although the weak-lensing mass threshold is, as expected, somewhat smaller for virialized lenses than for dark lenses 
(see Figure 7) their normalized distributions are not very different. Nonetheless, the sky density of virialized lenses is 10 — 80 
degree"^ and hence a factor of 4 — 10 larger than the sky density of dark lenses. This is because by redshifts oi z ~ 0.5 (where 
the distributions peak) a majority of objects in the mass range that can lens will have already virialized. 

Having computed the redshift distribution and sky density of dark and virialized lenses we now determine their relative 
abundances. The fraction of weak lenses that are caused by dark, non-virialized objects as a function of redshift is shown 
in the left panel of Figure 8. Out to z « 1 the fraction is nearly constant with about 20% of all weak lenses arising from 
dark objects. We again emphasize that the predicted abundance of dark lenses relative to virialized lenses is significant not 
because dark lenses comprise the lower-mass end of the mass function; on the contrary, virialized lenses have a lower mass 
threshold than dark lenses as shown in Figure 7. Rather, it is significant because according to the STHC model, a substantial 
fraction of cluster-mass objects are sufficiently overdense to produce a detectable lensing signal despite not having reached the 
virialization overdensity. For z > 1 the abundance of weak lenses of all types (both virialized and dark) drops off significantly. 
This, again, is because (zs) ~ 1 and because the evolution of overdensities has not yet had enough time to produce sufiiciently 
large overdensities. This is also illustrated in the right panel of Figure 8 where we show the fraction of virialized objects that 
can lens as a function of redshift. For 0.2<z<0.5 a large fraction of virialized objects with M > 5 x 1O^^M0 can produce a 
detectable weak-lensing signal but, for the same reason as above, by z = 1 this fraction is nearly zero. 

Finally we would like to point out that given the above results for the weak-lensing mass distribution, it is not surprising 
that in their study of weak lensing by low-mass gala^xy groups, Hoekstra et al. (2001) could only (just barely) detect a weak- 
lensing signal by stacking 50 such groups together. Namely, the groups in Hoekstra et al.'s sample, which were at a mean 
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Figure 6. The number counts of dark and virialized lenses for a truncated isothermal sphere {dashed curve), an Hernquist profile {dotted 
curve), and an NFW profile {solid curve). The ordinate gives the sky density of lenses at redshifts less than z. The top curve for a given 
density profile corresponds to the sky density of virialized lenses and the bottom curve to the sky density of dark lenses. 



redshift of 0.3, had a mean overdensity of only ~ 75 and a mean mass of just ~ 4 X 10" M0, assuming an isothermal density 
profile and using their measured value of ~ 275 km s~^ for the lensing-inferred velocity dispersion. Therefore, as Figure 7 
suggests, individual groups from their sample were neither massive enough nor sufficiently overdense to produce a detectable 
weak-lensing signal. 



5.3 The efFect of increasing the image size on the lensing signal 

In the above calculations wc assume that the lensing images arc 5 arcmin in radius, roughly the size of lensing maps to date. 
However, if a lens is relatively nearby or has a large radial extent it is possible that a large fraction of the total lensing signal 
is missed. This effect might be especially troublesome for the detection of dark lenses, given that they are not yet virialized 
and hence have larger radii. Wc now address this issue by determining the extent to which increasing the image size alters 
the predicted abundance of dark lenses. 

Before moving on however, we note that while we examined the predicted distribution and sky density of weak lenses for 
a variety of profiles, there is good reason to regard the NFW profile as the most plausible. For virialized lenses this is clearly 
the case as N-body simulations show that the halo density profiles are well fit by the NFW form. Though it is difficult to 
be as certain in the case of dark, non- virialized, lenses (N-body simulation fits to profiles have so far only been for virialized 
systems), because most of the dark lenses are well past turnaround (1 + 5^^ ^50) and because it is unlikely that the STHC 
model perfectly describes the evolution of overdensities all the way to virialization, assuming an NFW profile for dark lenses 
is a fair approximation. Furthermore, since virialization is expected to occur from inside-out, the centers of dark lenses, where 
most of the lensing signal is coming from (as we show quantitatively below), are likely near virialization and hence well 
described by the NFW profile. For these reasons (and also to avoid overly cluttered figures), the rest of the figures in this 
paper show results only for the NFW profile. To obtain approximate results for the other profiles, simply scale by the relative 
abundances shown in Figure 6. 

In Figure 9 we plot the cumulative signal-to-noise ratio and the fraction of the total signal as a function of the angular 
distance from the lens center for a lens with an NFW profile at redshift z = 0.3 with mass M = 1O"M0. Although the 
fraction of the signal that comes from within 5 arcmin is ~ 90% for lenses of overdensity 1 + (5^f; = 200 and 1 + S^];^ = 70 the 
lens with overdensity 200 requires an image size of just ~ 2 arcmin to be detectable {S/N = 5) while the lens with overdensity 
70 requires ~ 10 arcmin to be detectable. In general we find that in order to detect nearly all dark lenses with S/N > 5 in 
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Figure 7. The predicted weak-lensing mass distribution at redshift z = 0.5. Shown is the conioving number density as a function of mass 
of dark lenses {thick lines) and virialized lenses {thin lines) of mass M in the interval dlnM. Plotted for comparison is the virialized 
mass distribution, i.e.,^^(Af) x M {thin, dash-dot line). Since the minimum mass overdensity needed to produce a detectable lens (see 
Figure 4) is lowest for the NFW profile {solid curves), such a profile predicts a smaller weak-lensing mass threshold as compared to 
the Hernquist profile {dotted curve) and the truncated isothermal sphere {dashed curve). The two diamonds on each dark lens mass 
distribution curve mark the mass at which the minimum overdensity needed to produce a detectable lens is 275 and 100. The sharp 
lower-mass cutoff in the dark lens mass distribution is a consequence of the heaviside step-function nature of /dark- 



a given field the image area must be at least ~ 7r(15')^ w 0.2 degree^, as shown in Figure 10. Larger image sizes will not 
significantly increase the number of dark lenses detected as very little signal comes from radii larger than 15 arcmin. Also 
note that although lensing geometry favors a lens midway between observer and source, this effect is somewhat countered by 
the fact that the closer a weak lens is to the observer, the closer the source galaxy images pass to the highly overdense lens 
center (i.e., the solid angle subtended by the lens is larger). As a result, if the image size is large enough to enclose a large 
portion of the lens core, the lensing signal will be most strong when the lens-observer distance is slightly smaller than the 
lens-source distance. This accounts for the shift, shown in Figure 10, of the minimum of 1 -|- S^j^{z) toward smaller redshift 
as the image size is increased. 



5.4 Estimating ag from the abundance of weak lenses 

The present-day abundance of rich, X-ray clusters has been used to constrain the value of ag, the amplitude of mass fluctuations 
in spheres of radius 8h~^ Mpc (Evrard 1989; Henry & Arnaud 1991; White, Efstathiou & Frenk 1993; Viana & Liddle 1996; 
Eke, Cole & Frenk 1996; Kitayama & Suto 1997, Viana & Liddle 1999). In Figure 11 we show the extent to which the 
measured abundance of weak lenses (from future cosmic-shear surveys, say) can further constrain as. Here we have plotted 
the sky density of weak lenses (both virialized and dark) at the 95% confidence limits of ag given by Viana & Liddle (1999). 
Since weak lenses are produced by only relatively rare objects, their abundance is very sensitive to the value of as, suggesting 
the usefulness of weak lenses in measuring the amplitude of mass fluctuations. 

Another benefit of using weak lenses to measure as is their broad redshift distribution. In particular, a systematic 
uncertainty in measuring as by measuring rich cluster abundances is the degeneracy between as and Qm that arises from the 
limited range in redshift in which rich clusters are observed. To break this degeneracy substantial effort is made to measure 
not only the present-day rich cluster abundance but also the rich cluster abundance at higher redshifts {z m 0.3; e.g., Henry 
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Figure 8. Left panel: The fraction of weak lenses that are dark lenses as a function of redshift for a truncated isothermal sphere {dashed 
curve), an Hernquist profile {dotted curve), and an NFW profile {solid curve). The fraction is relatively constant between redshifts z = 0.2 
and z = 1.0, beyond which the abundance of both dark and virialized lenses drops to zero. The coarseness of the curves for z > 1 is 
an artifact of numerical noise that is a result of this drop off in both abundances. Right panel: The fraction of virialized objects with 
M > 5 X IO^^Mq that are able to weak lens as a function of redshift for the same density profiles as above. 

1997). This, in turn, gives an estimate of the evolution of the cluster mass function and hence an estimate of Q.m- However, 
because they are faint, detecting high redshift (z > 0.3) rich clusters is difficult. Weak lenses, on the other hand, do not suffer 
from this limitation and in fact are expected to have a broad redshift distribution and be most abundant at z ~ 0.5 (see 
Figure 5). As a result, detecting weaJs lenses provides an excellent means of measuring the evolution of the mass function and 
hence measuring fim- By thereby breaking the degeneracy between as and Q,m, weak lensing surveys are also well-suited to 
constrain the power-spectrum amplitude, erg. 



6 DISCUSSION AND CONCLUSIONS 

In this paper, we have calculated the abundance of dark and virialized lenses. This was accomplished by first expressing the 
lensing signal strength as a function of the dark-matter overdensity and redshift. Having determined the overdensity required 
to produce a detectable weak-lensing signal we used the STHC model to calculate the differential abundance of overdensities 
as a function of position along their evolutionary cycle. Overdensities whose lensing signal yielded S/N > 5 were divided into 
two classes: those with 1 + S^^ < 1 + S^^ were dark lenses while those with 1 + 3^^ > 1 -|- were virialized lenses. 

The distinction between dark and virialized lenses was based on the former being at an unrelaxed, and hence earlier, 
stage in the overdensity evolutionary cycle. This distinction is not axbitraxy but rather is expected to result in observational 
features that definitively separate the two classes of lenses. For instance, since dark lenses will typically have overdensities of 
1 + 5^^ ^ 100 while virialized lenses have 1 + 5^^ ^ 300 (see Figure 4) , the projected surface density of a dark lens is smaller 
than that in a virialized lens by a factor of 3^^"^ ~ 2. The sky density of galaxies in a dark lens will therefore be about two 
times smaller than in a virialized lens. As it is difficult to detect a significant galaxy overdensity for even a virialized, lensing, 
cluster at redshifts oi z ^ 0.5, it will be all the more difficult to do so for a dark cluster. Another distinctive observational 
feature expected of dark lenses is a low X-ray luminosity as compared with virialized lenses, a consequence of the X-ray 
luminosity function's steep dependence on total virialized mass. This effect might also account for the low X-ray luminosities 
observed by Postman et al. (2001) in three high-rcdshift clusters; namely, these objects are in fact proto-clusters that have 
not yet completely virialized. 

Although we considered a variety of density profiles in our calculations of the predicted distribution and sky density 
of weak lenses, as we noted in Section 5.3, there is good reason to regard the NFW profile as the most plausible form for 
both virialized and dark lenses. Nonetheless, while N-body simulations show that virialized systems are well fit by the NFW 
form, testing whether non-virialized, cluster-mass, halos in N-body simulations are also well described by the NFW profile is 
a worthwhile investigation that has not yet been performed. That said, we have shown that the redshift distribution of dark 
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Figure 9. Upper panels: The sigiial-to-noise ratio within an angular radius 9 from the lens center as a function of 9 for an NFW profile. 
Lower panels: The fraction of the total lensing signal that comes from outside the angular radius 9. Approximately 90% of the lensing 
signal eomes from a region smaller than the lensing halo radius. All four panels correspond to an object at redshift z = 0.3 and mass 
M = IQI^Mq. The left-hand-side plots are for a nonlinear overdensity of 200 and the right-hand-side plots for a nonlinear overdensity 
of 70. 



and virialized lenses for all the considered profiles is fairly broad with an average around z = 0.5 and a FWHM of Az Ri 0.5. 

The sky density of dark lenses for the NFW profile was calculated to be ~ 20 degree"^ (and ^ 10 degree"^ for an isothermal 
sphere profile) and should therefore be readily detectable by upcoming cosmic shear surveys. For virialized lenses, we found 
a sky density of ~ 80 degree^'^ assuming an NFW profile (and ~ 50 degree"^ for an isothermal sphere profile), a factor of 4 
to 5 larger than that of dark lenses. This difference is due to the fact that most of the weak lenses axe at redshift z w 0.5 and 
have masses of ~ 10^'' M© so that the majority are, according to the STHC model, virialized. It is important to note that 
while the aperture mass weight function used here was chosen to match a specific density profile, as we show in an upcoming 
paper (Weinberg & Kamionkowski 2002), a more general, non-optimal weight function, such as that given by Schneider et 
al. (1998), would lower the overall abundance of both types of lenses equally. The principle result of this paper, namely the 
expectation that ~ 10 — 20% of weaJj-lenses are dark, would not change. 

We find it encouraging that given the sky coverage of weak-lensing maps to date (~ 1000 arcmin^) and the average size 
of the individual lensing maps (~ 30 arcmin^ ), the number of dark lenses wo would expect to have seen is of order unity and 
thus consistent with the detection (Erben ct al. 2000, Miralles et al. 2002) of one or two dark lenses. Furthermore, in mock 
observations of numerical simulations. White, van Waerbeke, & Mackey (2001) showed that a weak-lensing search for clusters 
will likely suffer from serious line-of-sight projection effects due to the fact that clusters preferentially live in larger structures. 
These structures on larger scales, which perhaps correspond to 2a - 3a peaks in the primordial distribution, may well be the 
type of systems that wo find give rise to dark lenses. 

Finally, we have also shown that measuring the abundance of weak lenses can substantially help to constrain as, the 
rms mass fluctuation in spheres of radius 8h~^ Mpc. This is a consequence of the broad redshift distribution of weak lenses 
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Figure 10. The effect of increasing the image size on the abundance estimates. Plotted are the minimum overdensity needed to produce 
a detectable lens {left panel) and the sky density of weak lenses (right panel) as functions of redshift for image sizes of 0.02 degree^ 
{line-dot curve), 0.2 degree'^ {dashed curve), and 1 degree'^ {solid curve). The top three curves in the right panel correspond to virialized 
lenses and the bottom three to dark lenses. An NFW profile is assumed. Note that increasing the image size beyond 1 degree^ barely 
increases the predicted sky density since the signal becomes increasingly small at larger angular distances from the lens center (see Figure 
9). 

and the fact that they correspond to high-density peaks in the Gaussian primordial distribution. Cosmic-shear surveys, with 
their ability to detect cluster mass weak lenses over large areas of sky, should therefore provide a powerful new technique for 
determining the power-spectrum amplitude. 
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APPENDIX A: DERIVATION OF THE SIGNAL-TO-NOISE RELATION FOR VARIOUS DENSITY 
PROFILES 

Starting from equation (14), we derive the signal-to- noise relation for a point mass, a uniform-density sphere, a truncated 
isothermal sphere, an NFW profile and a Hernquist profile. Since these profiles are all axially-symmetric, {n){6) = k{6). 

(i) Point Mass: The dimcrisiorilcss mean surface mass density within a circle of radius 6 for a defiecting lens of point mass 
M at angular diameter distance Dd is 

where P = 6Dd. The quantity k, the dimcnsionlcss mean surface mass density on a circle of radius 9, is oc S{6), the Dirac 
delta function. Therefore, by equation (14), the signal-to- noise relation for a point mass is given by 

S_ V2{Z)M heoutY 1 (A2) 

N ae^crit^nDl (9out V V ft„ 7 ' ^ ' 

This can be expressed as a minimum mass needed to produce a detectable weak-lensing signal, which in useful units is 

n \-V2 
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where Dh ~ c/Hq is the Hubble distance. 

(ii) Uniform Density Sphere: Repeating the same procedure as above but for a sphere of uniform density p(r) = pc and 
mass M yields the following for the surface mass density (where we use the Abel integral equation to relate volume mass 
density to surface mass density); 



— / dz p{r) = ^^^/R^-P\ 

»itoo J -00 ^critoo 

-l^^\(P')2.P'dP' 



4pc ( P:^ - -{R^ - p2)3/2 



(A4) 



where dz is along the line of sight and R = {3M/4:-Kpc)^^^ is the radius of the sphere. The signal-to-noise ratio is then computed 
by solving equation (14) with the above relations for k and k. 

(iii) Truncated Isothermal Sphere: The radial density profile of am isothermal sphere is 

where cr„ is the line-of-sight velocity dispersion of the particles (i.e., galaxies) in the system. The surface mass density is then 

given by 

^=r=^4^ (A6) 

where 9e = cr^/GDdEcritoo- For a truncated isothermal sphere of mass M and radius R, M = drp{r)47rr^ = 2a^R/G, so 
that 

2R DcjScritoo ^ ^ 

Equation (14) then gives 

(iv) NFW Profile: The NFW density profile is given by 
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where Vb and ps are the scale radius and density, respectively. The mass within radius r is then 

M(r) = 47rp,rf (^ln{l + r/rs) - J^^jy^r) • (AlO) 

Bartelmann (1996) (see also Wright & Brainerd 2000) showed that the radial dependence of the tangential shear for an NFW 
profile is 

7nfw(a;) = k{x) - k{x) = ^"^^ g{x), (All) 
where x = ODa/ra and 
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The signal-to-noise ratio is then 
S_ _ 2V7m(Z) 

N ~ CTeEcritoo-Dd 



^out 

da;a;5(a;)2. (A13) 



There are thus three unknowns if given an overdensity of mass M: Vs, ps and R. We therefore need a third relation in 

addition to Equations (AID) and (A13) in order to break the degeneracy. It is obtained via the following conservation of 
energy argument, first put forth by Dalcanton et al. (1997) for the case of disk formation. 

Assume the mass profile before collapse is a uniform sphere of radius Ri and assume that at this initial stage the system's 
energy is entirely gravitational {E = —3GM^/5Ri). As noted by Dalcanton et al. (1997), this assumption is well motivated 
in the context of disk formation by the observed similarity between disk angular momentum distributions and the angular 
momentum distribution of a uniformly rotating sphere. It is natural to assume a similar initial condition occurs for systems 
at larger scales, i.e., cluster masses. As the overdensity collapses and approaches virializatiori, the mass distribution evolves 
into an NFW profile, as suggested by numerical simulations. At this stage the systems potential energy within a radius r is 

HV) = -Sn'Gpir! (l - ^^^) , (A14) 

where y = l + r/rs- Assuming the energy of the overdensity within Ri is conserved during collapse and that the system is near 
virialization so that E ~ \^\ /2 then gives Ri = 8.74rs. Since the truncation radius is given by the radius that contains mass 
M, by conservation of mass R = Ri = 8.74rs (i.e., though the mass is redistributed as the overdensity evolves the size of the 
sphere containing mass M is constant in time). When we include the effects of the cosmological constant in the conservation of 
energy argument there is little change in the result. The above relation between R and rs thus provides the sought after third 
equation needed to break the degeneracy between r^, ps and R. In an upcoming paper (Weinberg & Kamionkowski 2002) 
we show that the above approach yields concentration parameters that are slightly different from those obtained by N-body 
simulations (i.e.. Bullock et aJ. 2001). Nonetheless, the concentration parameters obtained by the two approaches predict a 
similar abundance of virialized lenses. Note that since the N-body simulations fit the concentration parameters to virialized 
objects, the above analytic approach must be used in order to compute the abundances of dark lenses, 
(v) Hemquist Profile: The Hernquist profile is given by 

P(^) = T-, — 77 (Alo) 

2n {r/rs){r + rsY 

where Vs is the scale radius and Moo is the mass enclosed at infinity. The mass within radius r is then 

2 

1 +r/' 

Using the Abel integral equation it can be shown that the dimensionless surface mass density for the Hernquist profile is 



M(r) = Moo I -^-^ 1 • (A16) 
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where x = ODa/rs and 
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The dimensionlcss surface mass density within x is then, 
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The signal-to-noise ratio is then obtained by inserting the above relations into ecjuation (14). 

As in the case of the NFW profile, given an overdensity of mass M, there are three unknowns. We therefore apply the same 
energy conservation argument as above, assuming the overdensity is initially a uniform density sphere of radius Ri and upon 
collapse relaxes to a Hernquist profile. The potential energy upon collapse is 



Hv) = - 



GMj 
6rs 



6t/" - 8y + 3 



(A20) 



where y = 1 -\- r jrs- Assuming energy conservation and a nearly virialized overdensity yields R = Ri = 3.2rs, allowing us to 
solve the signal-to-noise relation. 
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